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Abstract of EP01 96259 

Plant permitting a changing chemical medium to 
be controlled and regulated to arrive at an 
optimum value of a predetermined function of 
physical or chemical parameters representing the 
medium, the function Y being related to the 
values of the parameters by a relationship of the 
form: the said plant comprising: - sensors (5, 6) 
intended to supply the instantaneous value of 
each of the said parameters, - a scanning unit 
(26) allowing the value of the group of 
parameters to be recorded at uniform intervals, - 
a calculating unit (10) making it possible to 
determine the modification to be made to the said 
parameters to give an optimum value to the 
function Y, and - means for adjusting the said 
parameters controlled by the said calculating unit 
(10) to bring the value of the said parameters to 
the optimum. 
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(g) Procede de regulation de milieu chimique. 

<§) Installation permettant de contrfiler un milieu chimique 
evolutif et de le reguler, pour parvenir a une valeur optimale 
d'une fonction predeterminee de parametres physiques ou 
chimiques representatifs du milieu, la fonction Y etant reliee 
aux valeurs des parametres par une relation de la forme: 

Y = b 0 + b,X, + b 2 X 2 ... + b k X k + ... + b n X n + b, 2 X,X 2 + 
+ b k . iM X k .,X k + ...+ b kk X k * +...+ b nn X n * +... 

ladite installation comprenant: 

— des capteurs (5, 6) destines a fournir la valeur instanta- 
nee de chacun desdits parametres, 

— un organe de scrutation (46) permettant de relever perio- 
diquement {a valeur de I'ensemble des parametres, 

— un organe de calcul (10) permettant de determiner la mo- 
dification a apporter auxdits parametres pour donner une 
valeur optimale a la fonction Y, et 

— des moyens de reglage desdits parametres commandes 
par ledit organe de calcul (10) pour amener la valeur desdits 
parametres a I'optimum. 
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Procede et installation permettant d'une part, de connai- 
tre d'av ance a un instant donne le resultat auauel conduit 
necessai rement un milidu chxVniaue evb-Xuti^f^onn^ et . d' au- 
tre par t, d' assurer la regulation de ce milieu pour parve- 
nir a un resultat determine- fixe a l'avance. 

L* invention a pour objet un procede permettant, 
d'une part, de connaitre d'avance a un instant donne le 
resultat auquel conduit necessairement un milieu chimique 
evolutif donne et , d' autre part, d' assurer la regulation 
de ce milieu pour parvenir a un resultat determine fixe a 
l'avance. 

On rencontre dans l'industrie de nombreux milieux 
chimiques qui sont evolutif s de par la variation des para- 
metres qui les caracterisent , a savoir concentration en 
les divers cons tituants , pH . temperature et autres. 

A titre d'exemples de tels milieux, on peut citer 

- les bains acides de conversion chimique, 

- les bains alcalins ou acides de degraissage des metaux, 

- l'eau contenue dans des chaudieres et installations de 
chauf f age , 

et autres. 

L ' utilisation de ces milieux conduit a des resul- 
tats souhaites dans certains cas et non souhaites dans 
d a autres et qui peuvent se traduire dans le cas des exem- 
ples susdits respectivement : 

par des proprietes plus ou moins favorables des 
revetements obtenues par conversion chimique, 

- par un etat de proprete plus ou moins pousse des 
surfaces degraissees , 

; - par des etats d'entartrage variables des chau- 

dieres ou installations de chauf f age . 

Or, il importe a l'industriel de savoir a tout 
instant si le resultat obtenu a 1' instant en question par 
1 ' utilisation du milieu cor respondant est bien celui qui 
est attendu . 
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Mais Xa reponse est souvent longue et difficul- 
tueuse a obtenir. Ainsi, par exemple dans le cas d'un bain 
do conversion chimique, 

il faut disposer de 5 ou 6 semaines pour savoir 
si Xa resistance au brouiXXard saXin du revetement obtenu 
est superieure a 1000 heures ou non, 

- iX faut disposer d'un materiel couteux et encom- 
brant (microscope electronique a baXayage) pour connaitre 
Xa structure cristaiXine du revetement obtenu. 

Et ii est bien certain qu'une fois ces tests ef- 
fectues. c'est-a-dire par exemple 6 semaines pXus tard . Xe 
bain a pu evoXuer de maniere importante, sa composition 
pouvant aXors etre tres differente de ceXXe qu'ii avait au 
moment ou avait ete obtenu Xe revetement teste. IX s' en- 
suit que Xa modification qu'on aurait eventueXXement du 
apporter au bain a X* instant en question pour faire en 
sorte que Xe resultat obtenu soit ceXui souhaite, se trou- 
ve en fait depassee au moment ou Xe resuXtat est connu a 
X" issue des tests en question. 

IX existe deja des procedes de regulation de mi- 
lieux chimiques evolutifs. Mais, dans ces procedes, ii est 
d" usage de mesurer de fapon repetitive au moins l'un des 
parametres d'un teX milieu r de comparer chaque fois Xa va- 
Xeur mesuree a une vaXeur de consigne et d'agir sur Xe 
parametre en question de maniere a rapprocher Xa vaXeur 
mesuree suivante de Xa vaXeur de consigne, par exempXe en 
ramenant, si Xe parametre mesure est un constituant du 
miXieu evoiutif. Xa concentration en ce constituant a Xa 
vaXeur de consigne par apport de produits de regeneration 
appropries . 

ii est bien certain que, dans un milieu chimi- 
que evoiutif du genre en question, ces operations ne peu- 
vent garantir la Constance des resultats que I'on souhaite 
obtenir. 

Le probleme que la Societe Demanderesse s'est done 
propose de resoudre est celui de mettre au point un precede 
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propre a permettre a 1 1 utilisateur : 

d'une part, de connaitre, a partir de la mesure 
de divers parametres evolutif s d'un milieu chimique evo- 
lutif, les particularites du resultat auquel conduit le 
milieu en question a 1' instant choisi, 

d' autre part, au cas ou la valeur trouvee par 
cette mesure montre que le resultat en question ne corres- 
pond pas a celui souhaite, de determiner tout d'abord un 
ensemble de modifications des variables ou parametres du 
milieu. propre a faire en sorte que le resultat auquel 
conduira le milieu ainsi modifie corresponde a celui sou- 
haite, puis apporter cet ensemble de modifications au 
milieu chimique evolutif en question, par action sur les 
parametres concernes . 

Et elle a eu le merite d'elaborer un procede per- 
mettant d'atteindre ce but, ce procede etant caracterise 
par le fait 

que, dans le cas d'un milieu chimique evolutif 
donne et pour au moins l'un des resultats auxquels conduit 
ce milieu, on etablit une equation dite pre visionnelle du 
type 

Y - b Q ♦ b lXl ♦ b 2 X 2 ... . b k X k ..... b n X n ♦ b 12 X lXz .... 

+ b k-1.k X k-1*k b kk X k 2 b nn X n 2 

dans laquelle 

- Y represente le resultat auquel conduit le milieu. 

- X« . X o .... X. , . . .X representent les variables ou 

it k n 

parametres du milieu chimique evolutif dont depend le 

resultat Y et 

- b_ . b. . . .b , . . . , b sont des coefficients par 
u I k nn 

rapport auxquels s'exprime la linearite de 1" equation, la 
valeur de ces coefficients etant determinee par des cal- 
culs reposant sur des methodes de regression en elles- 
memes connues , 

que l'on mesure de maniere repetitive chacun 
desdits variables ou parametres par autant de sondes ou 
capteurs et 
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que. d'une part, a un instant t donne, c'est-a- 
dire a 1" issue d'une des mesures repetitives, on caXcule 
la valeur de chacune des equations Y en utilisant les va- 
leurs mesurees pour les differents variables a 1* instant t, 
5 - que, d' autre part, e ven tuellement et au cas ou 

la valeur trouvee pour I'une ou l'autre des equations Y ne 
correspond pas au resultat souhaite, on determine par le 
calcul un ensemble de variations ou modifications a impar- 
tir aux . variables ou parametres, cet ensemble de varia- 

10 tions donnant pour Y la valeur souhaitee et constituant, 
parmi toutes les solutions possibles de 1' equation previ- 
sionnelle Y, celle qui est la plus economique et la plus 
-Facile, a realiser r et 

que r d'autre part encore, on agit individ uelle- 

T5 mervt sur autant d'organes de reglage que de parametres, 
fournissant a chacun de ces organes de reglage un signal 
propre a amener la variable ou parametre qui lui est asso- 
cie a la valeur cor respondant a la susdite solution de 
1* equation previsionnelle Y la plus economique et la plus 

20 -Facile a obtenir. 

Dans un mode de realisation avantageux, on effec- 
tue les calculs susdits a l'aide d'une calculatrice nume- 
rique que 1'on munit de bases de donnees appropriees . 

L ' invention vise encore d'autres dispositions qui 

25 s'utilisent de preference en mime temps et dont il sera 
plus explicitement question ci-apres. 

Et elle pourra , de toute facon, etre bien comprise 
a l'aide du complement de description, du dessin y relatif 
et des exemples , lesdits complement de description et 

30 exemples concernant des modes de realisation avantageux. 

La figure 1 de ces dessins est un diagramme rela- 
tif a une equation Y et la figure 2 est une representation 
schematique d'une installation propre a la mise en oeuvre 
du procede selon 1' invention. 

35 Les figures 3a et 3b prises dans leur ensemble, 

d'une part. et 4a, 4b et 4c prises dans leur ensemble, 
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d' autre part, representent chacune un logigramme relatif a 
la mise en oeuvre du procede conforms a l' invention. 

La figure 5 est un schema synoptique d'un exemple 
d ' installation permettant de mettre en oeuvre 1" invention. 

Se proposant, par consequent, de se donner les mo- 
yens permettant, d'une part, de connaitre a l'avance, a un 
instant donne, le resultat auquel conduit necessairement 
un milieu chimique evolutif et , d' autre part, d' assurer la 
regulation de ce milieu pour parvenir a un resultat deter- 
mine fixe a l'avance, on s'y prend comme suit ou de fagon 
equivalente . 

Pour un milieu chimique evolutif donne. il peut 
etre necessaire de connaitre a tout instant le resultat ou 
propriete Y ^ auquel il conduit, ce resultat ou propriete 
etant la fonction d'un certain n ombre de variables X 1 , X 2# 

X^* du milieu et dont depend le resultat, Y et les 

variables X^ . X^ , ... X i etant lies par 1' equation previ- 
sionnelle 

Y« 55 *> n + b 4 X, ♦ b.X. b. . X.X. ♦ b. .X. 2 ... 

1 0 11 XX 1} X □ x.xx 

dans laquelle b Q , b 1 ... b^ ... sont des coefficients, au 
depart inconnus et determines de la maniere qui sera indi- 
quee plus loin, la valeur Y^ effective traduisant a un mo- 
ment donne le resultat auquel conduit le milieu pris dans 
sa composition correspondant a ce meme moment. 

En admettant que le resultat en question doit etre 
par exemple egal ou superieur a la valeur N, il faut que. 
pour que le milieu a 1' instant envisage puisse etre consi- 
dere comme convenable, qu'a cet instant la mesure des va- 
riables X 4J , X 2 ... X i fournisse pour une valeur 

Y > N 

la mesure des variables X^ , ... X^ pouvant etre effec- 

tuee par autant de capteurs ou sondes plongeant dans le 
milieu et affectes respectivement a chacune des variables 
X ^ , ^2 • • • X ^ • 

Si Y. ainsi mesuree est * N, le milieu conduit au 
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resultat souhaite ; sinon il faut apporter des modifica- 
tions au milieu. 

Mais avant toute chose, il faut bien entendu eta- 
blir l'equation Y et, par consequent, 

- selectionner les variables X. que l'on considers comme 
etant caracter istiques du milieu (il peut s'agir de la 
concentration en certains constituaots, de la temperature, 
du pH et autres ) r 

- determiner les valeurs des coefficients b Q . . . b^ .... 



Le choix des variables est determinant pour 1'ob- 
tention d'une equation previsionnelle de bonne qualite. 
1'oubli d'un parametre influent conduisant a l'etablisse- 
ment d'une equation erronee. Pour chaque probleme pose, le 
choix des parametres resulte des connais sances anterieures 
15 cfe l'homme de l'art. 

La determination des coefficients b^^ necessite le 
recours a des calculs reposant sur des methodes de regres- 
sion classiques (sur lesquelles il ne sera pas insists ici 
mais dont le principe peut etre trouve dans l'ouvrage de 
20 L. LEBART, A. MORINEAU, i.P. FENELON, "Traitement des 
cfonnees statistiques t methodes et programmes", DUNOD ) en 
utilisant les donnees fournies par un certain nombre d" ex- 
periences. Du point de vue pratique, on precede pour ce 
faire de la maniere suivante : 
25 - on choisit la ou les proprietes a me surer , 

oh choisit les parametres experimentaux ou va- 
riables qui vont decrire l'equation, 

- on choisit pour chaque parametre un domaine de 
variations qui est fonction du probleme pose et de la con- 

30 naissance anterieure de 1 * experimentation , 

- on choisit un modele d ' equation simple, le pre- 
mier modele choisi pouvant etre du type lineaire : 

Y = b + b X + + b. X. ♦ + b X 

10 11 ii n n 



ou lineaire avec interactions 

Y 1 ~ b n * b i X i b * x -i + b - x *••♦ b x + -- + *>. -X.X.+ , 

1 0 11 il 3 3 nn 13 13 
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- on realise le nombre d ' experiences defini par le 
choix du modele et on releve pour chacune d'elles la va- 
leur de la reponse Y obtenue, 

- on calcule les coefficients . . . b^ • - - du modele. 
Le modele est alors teste. Ce modele etant choisi 

"a priori", il est en effet indispensable de le tester 
pour 1* accepter et l'utiliser. ou bien au contraire le 
refuser. Pour cela , on realise quelques experiences sup- 
plementaires pour lesquelles les valeurs de chacun des 
parametres sont fixees simult anement a egale distance des 
extremites de leur domaine de variations et la reponse 
obtenue est comparee a la reponse calculee par le modele. 
Si les deux resultats sont voisins, le modele est accepte. 
Dans le cas contraire, on complique le modele par adjonc- 
tion de termes supplementaires, par exemple des termes 
quadratiques . 

L* equation pr evi sionnelle representative du nou- 
veau modele s'ecrit alors : 

Y i ■ b o * b i x i + - + b x x x + *- + b ia x iV--- 
+ b ii x i 2 * b jj x j 2 * ••• 

On realise le supplement d ' experiences necessaires 
au calcul des coefficients des termes quadratiques et le 
modele est de nouveau teste. 

Ainsi, de proche en proche, par une demarche se- 
quentielle, on affine le modele jusqu'a ce que la previ- 
sion soit de bonne qualite. 

L * equation ainsi etablie est representee par une 
hypersurface dans un espace de n+1 dimensions pour n fac- 
teur s . 

Pour materialiser 1' incidence des divers facteurs. 
on projette 1 * hyper surf ace successivement sur des sous- 
espaces a trois dimensions dont deux sont constitutes par 
deux variables que 1 * on desire etudier et la troisieme par 
la reponse Y, les variables restantes etant fixees a des 
valeurs constantes ; ceci se traduit par des courbes de 



0196259 



niveau dans le pXan des deux variables, la reponse etant 
materialisee par les valeurs des courbes de niveau ; ainsi 
on peut erivisager de projeter 1 ' hypersurf ace representa- 
tive de I 4 equation Y ^ sur le plan des parametres X 3 et X^; 

5 a titre d'exemple, on peut imaginer qu'Y represente la 

i 

resistance au brouillard salin d'une surface metallique 
qui a ete soumise au prealable a un traitement de conver- 
sion chimique et a un revetement par peinture, X 3 repre- 
sentant par exemple la teneur du milieu evolutif , c*est-a- 

10 dire du bain de conversion, en ions Zn + * et X la teneur 

3- 

du meme bain en ions PO , , les variables restantes X.. et 

X^ etant fixees a des valeurs constantes. 

Dans ce cas r l'hypersurface se traduit dans le 

plan de-fini par X_ , X, par un certain nombre de courbes 

j w 

15 dites "courbes isorepoases" ou courbes de niveau et desi- 
gnees par C^ t . ..C^, cfiacune de ces courbes correspon- 
dent a une valeur donnee de la resistance au brouillard 
salin et chaque point de chacune des courbes representant 
un couple de valeurs de X 3 et X^ qui conduit a une surface 

20 metallique ayant la resistance au brouillard salin corres- 
pondant a la courbe. 

A la -figure T , on a represente un diagramme (pro- 
jection de l'hypersurface sur le plan X^. X^ > montrant les 
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courbes isoreponses suivantes i 
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Si on admet que, pour etre satisfaisante, la re- 
sistance au brouillard salin des surfaces traitees doit 
35 etre de 950 a 1000 h, tous les couples de valeurs de X 3 et 
X^ correspondant a des points situes sur l'une des courbes 
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C- et C Q ou dans la surface comprise entre ces courbes 

conduisent au resultat recherche. 

En d'autres termes si, a un moment donne, la mesu- 

re de X- et X. conduit a un point P repondant a cette 
3 4 i 

5 condition. la resistance au brouillard salin de la piece 
traitee sera satisfaisante ; si. par contre, elle conduit 
a un point P 0 ne repondant pas a cette condition, c'est-a- 
dire situe en dehors du domaine defini par les courbes C ? 
et C a , le resultat ne sera pas satisfaisant et il s'agira 

8 

10 d'apporter au milieu evolutif , c'est-a-dire au bain de 

conversion, une modification des concentrations X 3 et X^ 

qui ramene a 1'interieur dudit domaine le point situe en* 

, hors de celui-ci. 
de 

Et la modification retenue parmi toutes les modi- 
15 fications possibles conduisant a ce resultat devra etre 
celle qui soit la plus facile et la plus economique a rea- 
liser . 

Cette determination est realisee a chaque instant 

en effectuant les calculs necessaires pour determiner pour 

20 chacune des differentes variables concernees et non pas 

seulement X^ et X, , de chacune des equations previsionnel- 
3 4 

les . Y rt . . . Y la modification la plus simple et la 

1 2 n 

plus economique parmi toutes les modifications possibles 
conduisant a un resultat satisfaisant. 
25 Par exemple, on peut effectuer ces calculs a 

l'aide d * une calculatrice munie des bases de donnees ne- 
cessaires, notamment les donnees relatives aux equations 

previsionnelles Y 4 , Y 0 . . . Y . . , . 

1 2 n 

Pour pouvoir effectuer 1' ensemble des calculs, il 
30 convient de disposer a 1* instant donne pour chaque equa- 
tion previsionnelle de 1* ensemble des valeurs des parame- 
tres correspondants , ce qui peut etre realise dans le 
cadre des mesures repetitives dont il a ete question plus 
haut par des capteurs ou sondes plongeant dans le milieu 
35 chimique evolutif considere. 

La modification est alors realisee par 1' envoi a 
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un ensemble de dispositif s aptes a agir sur les divers 
parametres, de signaux propres a dedencher un ensemble de 
reponses de la part de ces dispositifs, reponses qui se 
traduisent au niveau de chaque parametre par un changement 
5 precis, 1" ensemble de ces changements conduisant a la sus- 
dite modification, selectionnee et determinee par les cal- 
culs . 

Pour la commodite de l\expose, on admet dans ce 
qui suit et en particulier dans les examples que le milieu 

to chimique evolutif pris en consideration est une sequence 
de traitements de phosphatation de bain de conversion chi- 
mique au zinc, etant bien entendu que les principes qui 
vont etre developpes en rapport avec cet exemple particu- 
lier de milieu chimique evolutif se transposent aisement a 

15 tout autre milieu chimique evolutif en tenant compte des 
considerations plus generales qui precedent. 

II est connu que les bains de conversion en ques- 
tion comprennent des ions nickel, zinc, phosphate et sont 
acceleres par des nitrites et/ou des chlorates. 

20 Par ailleurs, la sequence de traitements peut com- 

prendre de fapon classique, par exemple, un traxtemertt de 
conditionnement prealable de la surface et un traitement 
subsequent consistant en un rinpage passivant r il va de 
soi que la sequence peut comprendre d'autres traitements, 

25 mais ceux-ci ne sont pas pris en consideration ici pour la 
commodite de la description. 

La figure 2 montre schema tiquement une installa- 
tion propre a mettre en oeuvre le procede conforme a 1" in- 
vention . 

30 Cette installation comprend 

trois enceintes representees schematiquement en 
1 . 2 et 3 et renf ermant respectivement les bains de con- 
ditionnement, de conversion proprement dite et de rinpage 
passivant , 

35 - six sondes ou capteurs representes schematique- 

ment en 4 a 9, constitutes par des pH-metres, des electro- 
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des selectives. des thermometres et autres et plongeant 
respectivement , pour ce qui est de la sonde 4, dans 1* en- 
ceinte 1, pour ce qui est des sondes 5, 6 et 7 , dans 
l'enceinte 2 et, pour ce qui est des sondes 8 et 9, dans 
1 ' enceinte 3 ,. 

une calculatrice 10 dans laquelle ont ete d'une 
part introduites sous forme de programme les algorithmes 
de resolution des diverses equations Y n envisagees, d' au- 
tre part memorisees les donnees permanentes et coeffi- 
cients, cette calculatrice etant reliee aux sondes 4 a 9 
dont elle peut par consequent recevoir les valeurs mesu- 
rees par 1 * intermediaire d' interfaces comportant des con- 
vertisseurs A/N. 

une chaine de commande 1 1 reliee a la calcula- 
trice 10 et propre a commander, sous 1" influence de si- 
gnaux emis par la calculatrice et par 1 ' intermediaire de 
dispositifs s chema tiquement representes par des fleches 
12, 13 et 14, 1 ' alimentation en produits notamment chimi- 
ques appropries des diverses enceintes 1 , 2 et 3 et 

eventuellement une imprimante 15 fournissant en 
clair les variations des diverses equations Y n sur les- 
quelles on entend jouer et dont les donnees ont ete four- 
nies a la calculatrice . 

En rapport avec la susdite calculatrice, on a in- 
dique sur les figures 3a-3c et 4a-4c des logigrammes mon- 
trant l'une des approches possibles pour, tout d'abord, 
charger en memoire de la calculatrice des courbes d'isore- 
ponses sous forme discrete, d" autre part, rechercher 1 ' op- 
timum de la reponse Y par la methode du gradient. 
Trace et memorisation des courbes isoreoonses 

Le cas illustre par le logigramme des figures 3a- 
3c correspond a celui d* equations polynomiales dont cer- 
tains termes sont des fonctions d'un des parametres et 
d' autres des fonctions de deux parametres. 

On commence par indiquer, en memoire de travail de 
la calculatrice, le nombre N de facteurs experiment aux 
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X(i) qui interviennent dans Xa determination de Y. On in- 
troduit ensuite I et J , qui sont Xes deux valeurs de i 
constituant chacun des couples pour XequeX on va tracer 
des courbes d'isoreponses. Au cours du derouXement du 
5 programme, I et J prendront evidemment tous Xes coupXes de 
vaXeurs possibXes entre 1 et N, a X'exception des coupXes 
pour XesqueXs I=J. et pour chaque coupXe (X(I), X(J)), on 
affectera, aux N-2 variables restantes, tous Xes jeux de 
vaXeurs constantes . 
10 Une . fois ce chargement effectue, on introduit 

(cadre 20) X* equation de Xa fonction reponse que r pour 
simpXi-Fier. on a suppose etre de Xa -Forme Xineaire avec 
interactions : 

i = N N 
15 Y = bo + E b 4 XCi) «- Z b_ X(i) XCj) + . 

i=1 
3*i 



20 



25 



30 



35 



i=1 1 i=1 13 



Avant de tracer Xes courbes, on va determiner Xa 
vaXeur de Y pour tous Xes points d*une griXXe rectanguXai- 
re correspondent a des vaXeurs de XCi) et X(j) reguliere- 
ment reparties. Pour que Xes griXXes correspondant a tous 
Xes coupXes de vaXeurs aient Xe meme nombre de points, iX 
faut effectuer une normalisation et un centrage tenant 
compte de Xa pXage de variations possibXes pour chacun des 
parametres X. CeXa revient a faire varier Xe nombre repre- 
sentant chaque parametre X entre -1 et + f f ce qui impXique 
d* introduire un facteur d'echelle. Pour cela, on introduit 
en memo ire Xes vaXeurs Ximites des parametres (cadre 22) 
puis on determine, pour chaque coupXe de parametres cor- 
respondant a une courbe d'isoreponse, Xe -Facteur d'echeX- 
Xe. Le cadre 24 correspond a une griXXe a P points pour 
XCI) et Q points pour X(J). correspondant a un premier jeu 
de vaXeurs constantes de tous Xes X autres que XCI) et 
XC J) . 

Le trace des courbes peut aXors etre effectue a 
partir des vaXeurs minimaXes prevues pour X(I) et X(J) 
(par exempXe temperature ambiante, d'une part, concentra- 
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tion nulle, d'autre part). Apres initialisation a ces 
valeurs en 26, on effectue le calcul de tous les points en 
sequence, en mettant en oeuvre deux boucles imbriquees. Le 
cadre montre en 28 correspond a la determination de tous 
les points de la grille correspondant a une valeur donnee 
d'un des parametres avant de passer a une nouvelle valeur 
de ce parametre pour recommencer le calcul. 

On determine ainsi les valeurs prises par Y pour 
tous les points de la grille. On trace ensuite les courbes 
d ' isoreponses par une interpolation dont les parametres 
sont introduits en 30- Le pas d'ecartement des isoreponses 
(valeurs de Y pour lesquelles on veut disposer des cour- 
bes) sera choisi en fonction de la finesse d'ajustement 
que l'on recherche. 

L* operation sera effectuee a partir des points de 
la grille, en utilisant un programme d ' interpolation qui 
peut etre classique. 

A I'issue de 1'operation, schematisee en 32, tou- 
tes les valeurs d ' isoreponses seront affichees et memori- 
sees . 

La meme operation sera repetee pour tous les cou- 
ples X(I), XI J) avec tous les jeux de valeurs constantes 
affectees aux parametres fixes. 

Une fois les donnees ainsi stockees dans la calcu- 
latrice, cette derniere peut etre utilisee pour optimiser 
la fonction Y par la methode du gradient. 
Optimisation de Y oar la met hode du gradient 

Le logigramme montre en figures 4a a 4c reprend 
une fraction du logigramme precedent. Mais il ne met pas 
en oeuvre cette partie commune pour la totalite du domaine 
de variation des parametres X, mais seulement autour du 
point experimental, au lieu de repartir des valeurs mini- 
males . 

Une installation fonctionnant suivant le logigram- 
me des figures 4a a 4c constitue done un ensemble autonome. 
Ce logigramme prevoit des contraintes, visant a diminuer 
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la duree des calculs et a eviter des reglages trop fre- 
quents des parametres. Le parcours des boucles de calcul 
est arrete lorsque le gain sur Y obtenu par un parcours 
supplements ire est inferieur a une valeur predeterminee. 
5 De plus, les valeurs des reponses doivent rester a l'inte- 
rieur d' un domaine de validite defini par un rayon R. 

La -Figure 4a montre les operations d ' initialisa- 
tion du programme. On introduit H et l'equation de la 
fonction-reponse Y, c'est-a-dire des donnees identiques a 

to celles utilisees au debut du programme de trace de courbes. 
Pour appliquer la metfiode du gradient, on introduit ensuite 
les deriyees partielles D[X( i )] = &YfbX ± de Y par rapport 
a un premier parametre X^^ (cadre 34 > puis les valeurs li- 
mites et l*algorithme du gradient (cadre 36) avec une va- 

15 leur predeterminee, la meme pour tous les X. On determine 
ensuite en 38 (figure 4b) le sens de 1' influence sur Y 
d'une premiere variation de 1" ensemble des parametres 
avant d ' ef f ectuer la recherche de 1* optimum, par applica- 
tion de 1'algorithme du gradient. II est important de 

20 remarquer que le logigramme reduit le coefficient multi- 
plicateur applique a la derivee lorsque la derniere itera- 
tion precedente amene a depasser le rayon de validite R et 
recommence le meme calcul, mais avec la valeur reduite. 
Dans l'exemple represents r cette reduction est faite dans 

25 un rapport 2. Si le nouveau calcul conduit encore a depas- 
ser le rayon R, une nouvelle reduction dans le meme rap- 
port est appliquee : la precision est ainsi tres accrue. 

Le logigramme linvite par ailleurs le nombre de 
parcours de la boucle en arretant les iterations des que 

30 l*on arrive a proximite de 1* optimum, ce qui se traduit 
par le fait que la difference entre deux resultats succes- 
sifs est inferieure a un seuil DY ou que, a la fois, on a 
depasse le rayon R et que le coefficient L applique a la 
derivee dans 1'algorithme du gradient est tombe au-dessous 

35 d'une valeur determinee (test 44, figure 4c). 
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La figure 5 montre, a titre de simple example , une 
installation susceptible de mettre en oeuvre le procede 
qui vient d'etre defini. Sur cette figure, les organes 
correspondant a ceux de la figure 2 sont designes par le 
mime numero de reference, 

Dans la realisation rnontree en figure 5. un cir- 
cuit 46 de scrutation periodique echantillonne , par exem- 
pie une fois par heure, les parametres X, qui dans le cas 
illustre sont au nombre de deux et fournis respecti vement 
par les sondes 5 et 6 . La valeur echantillonnee est appli- 
quee a un convertisseur analogique-numerique 48 et envoyee 
a un organe de calcul 10, qui peut etre un microproces seur 
ou un automate programmable. Les nouvelles valeur s de con- 
signe des parametres X elaborees par 1* organe de calcul 10 
sont envoyees a des cartes 50 et 52 de reglage des parame- 
tres. Dans les systemes a cinetique lente, se pretant mal 
a une regulation, la carte 50 pourra etre prevue pour ela- 
borer, a partir de donnees connues sur le systeme, de la 
valeur de consigne qui vient d'etre fournie par 1' organe 
de calcul 10 et de la valeur reelle fournie par la sonde 
5, la quantite de reactif a ajouter au systeme pour obte- 
nir la valeur de consigne. L* operation de reglage peut 
alors etre effectuee par commande d'ouverture d'une vanne 
52 pendant une duree determinee. Dans le cas ou la cineti- 
que est plus rapide, on peut utiliser une carte 52 qui 
repoit egalement, en vue d'une regulation, un signal de 
sortie provenant du convertisseur 48, le circuit de scru- 
tation 46 fournissant alors en permanence la valeur du 
parametre correspondant. afin de permettre d ' en suivre 
Involution. La carte 52 arrete alors l ouverture de la 
vanne correspondante 54 des que la valeur de consigne est 
atteinte. 

C'est en mettant en oeuvre les procede et instal- 
lation conformes a 1 ' invention qu * on a realise les exem- 
ples suivants. 
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gyEMPLE 1 - Determination des equations previsionnelles 

representatives de deux proprietes d'une sur- 
face metallique traitee dans un bain de con- 
version chimique. 
5 Dans cet exemple , il a ete choisi de modeliser et 

done de prevoir deux proprietes d'un revetement de conver- 
sion couramment mesurees en controle dans les traitements 
de phosphatation, a savoir : 

2 

- le poids de couche exprime en g/m , 
10 - Xa densite de structure cristaXXine exprimee en 

nombre de cristaux par unite de surface et mesuree a X* ai- 
de d r une photographie prise au microscope electronique a 
balayage avec un grossissement de 1500 fois . 

Cinq variabXes X. a X_ du bain de conversion ont 
t5 ete prises en consideration : 

: concentration en ions nickel 
X 2 : concentration en ions chlorate 
X 3 : concentration en ions zinc 
X^ : concentration en ions phosphate 
20 x 5 : valeur de X'acidite Xibre ou AcL 

etant entendu que , par AcL , on design© Xa vaXeur obtenue 
exprimee en ml de NaOH N/10 sur un echantillon de 10 mX du 
bain dose par NaOH N/10 au virage de X * heXianthine . 

Des eprouvettes en toXe d'acier lamine a froid de 
25 qualite courante sont traitees suivant une gamme au trempe 
a base de produits commercialises par Xa Societe Demande- 
resse et couramment employes dans 1* Industrie Automobile 
pour la preparation des carrosseries avant peinture. 

Les etapes et produits mis en oeuvre comprennent 
30 - un degraissage alcalin a chaud avec utilisation 

des agents degraissants a base de sels alcalins et de 
tensio-actifs commercialises par la Societe Demanderesse 
respectivement sous les marques de fabrique " R IDOL INE 1550 
CF/2" et "RIDOSOL 550 CF" , 
35 - un rinpage -froid a l'eau courante, 

un conditionnement de surface avec utilisation 
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d * un agent d'affinage a base de sel de titane commerciali- 
se par la Societe Demanderesse sous la marque de fabrique 
"FIXODINE 5", 

un traitement de conversion chimique avec uti- 
5 lisation d'une solution acide a base de phosphate de zinc 
et autres ions classiques commer cialisee par la Societe 
Oemanderesse sous la marque de fabrique "GRANODINE 908", 

- un rinpage a froid, 

- un rinpage passivant avec utilisation d'un agent 
10 a base d'ions chrome commercialise par la Societe Demande- 
resse sous la marque de fabrique " DEOXYL YTE 4 1 M , 

- un etuvage. 

Les cinq variables du bain de conversion chimique 
qui sont prises en consideration sont mesurees par tout 
15 moyen classique habituellement utilise : dosage chimique, 
pH-metrie, potentiometrie , conductivity ionique et autres. 

On effectue 16 experiences pour determiner les 
constantes de 1* equation previsionnelle correspondante qui 
s ' ecrit : 



20 



25 



30 



35 



Y 1 3 b o * b i 


X 1 


2 2 


+ 


b 3 X 3 = b 4 X 4 + b 5 X 


+ b 13 X 1 X 3 




b 14 X 1 X 4 


+ 


b 15 X 1 X 5 * b 23 X 2 X 3 


* b 25 X 2 X 5 


+ 


b 34 X 3 X 4 




b 35 X 3 X 5 + b 45 X 4 X 5 



b 12 X 1 X 2 
b 24 X 2 X 4 



et on mesure pour chaque experience 

- le poids de couche de chaque eprouvette traitee. 
par dissolution du revetement dans une solution d* acide 
chromique a chaud (a 10 Z, 60 # C, 15 minutes), 

le nombre de cristaux presents a l'interieur 
d'une fenetre de 7,5 cm de diametre de chaque eprouvette 
traitee, ces cristaux etant comptes a partir d'une photo- 
graphie prise au microscope electronique a balayage avec 
un grossissement de 1500 fois. 

Pour determiner les car acteristiques de ces expe- 
riences de fagon a pouvoir determiner les susdites cons- 
tantes, on procede comme suit : 

Pour chacun des cinq facteurs pris en considera- 
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tion , cm definit un domains de variations comme suit : 
0 . 6 g/1 S Zn ++ < 1.4 g/1 
15 g/1 * P0, 3 ~ * 23 g/1 
0.3 g/1 * ClOj- S 0.9 g/1 
5 0.5 g/1 * Ni + * 0.8 g/1 

0,9 ml NaOH N/10 S AcL *1,3 ml NaOH N/10, 

Compte tenu de 1' equation precedente, il -Faut de- 
terminer d (d»1 1 coefficients lies aux termes du 1er degre 
2 

et aux termes rectangles (p representant le nombre de pa- 
10 rametres avec, dans le cas present p = 5). auxquels vient 
s'ajouter le terme constant- Cela conduit a realiser les 
16 experiences figurant dans le tableau suivant pour les- 
quelles chacun des parametres ne prend que les valeurs 
extremes de son domaine de variation. 

15 



25 



35 



r- «, _ 

bxp . 

nb, 


NX 


cio 3 - 


£ n 




A — . 1 V 

ACL T ^ 


p . d • c . 
(a/m2) 


Y 

1 n . c . 


1 


0 . 5g/l 


0 . 3g/X 


0 ,8g/X 


1 5g/X 


1 , 3ml 


1 ,74 


1 47 


2 


0 . 5g/l 


0.9g/l 


O, 8g/X 


15g/X 


0 . 8ml 


1.14 


160 


3 


0 , 8g/l 


0.9g/X 


0. 8g/l 


15g/X 


1 . 3ml 


1 r 19 


1 68 


4 


0,5g/l 


o.3g/i 


1 ,4g/X 


15g/l 


0 , 8ml 


1,61 


1 17 


5 


0 , 5g/ 1 


0 . 9g/X 


1 . 4g/l 


15g/X 


1 , 3ml 


1 .78 


91 


6 


0,8 g/1 


0 . 3g/X 


0.8g/X 


23g/X 


1 . 3ml 


1 , 39 


128 


7 


0,5g/l 


0 .3g/X 


1 ,4g/X 


23g/X 


1 . 3ml 


1 . 86 


101 


8 


0 , 8g/l 


0.9g/l 


1 ,4g/X 


15g/l 


o , 8ml 


1 . 29 


81 


9 


0.8g/l 


0 . 9g/X 


0 . 8g/X 


23g/X 


0 , 8ml 


1 . 19 


205 


10 


0.8g/l 


0 .9g/X 


1 ,4g/l 


23g/X 


1 . 3ml 


1 .41 


66 


1 1 


G,5g/1 


0.9g/X 


1 ,*g/i 


23g/l 


0 . 8ml 


1 .61 


88 


1 2 


0 , 5g/l 


0,9g/X 


0.8g/X 


23g/l 


1 . 3ml 


1 .36 


157 


13 


0,5g/l 


0,3g/l 


0 . 8g/X 


23g/l 


0 . 8ml 


1 r 34 


126 


14 


0 . 8g/l 


0.3g/X 


0 , 8g/l 


15g/l 


0 , 8ml 


1 , 39 


161 


1 5 


0,8g/l 


0.3g/X 


1 ,*g/x 


23g/l 


0 , 8ml 


1.81 


97 


1 S 


0,8g/l 


0.3g/X 


0.8g/l 


15g/l 


0 . 8ml 


1.91 


90 
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Pour s'affranchir des unites de facteurs, pour 
simplifier les calculs et les interpretations et pour com- 
parer les effets des facteurs entre eux, les valeurs de 
chaque facteur sont centrees et reduites. 
5 Le calcul des coefficients b^ se fait alors par la 

methode des moindres carres decrite en detail dans l'ou- 
vrage cite precedemment . Seuls sont pris en compte dans 
1'equation, les coefficients les plus influents. 

Apres traitement des donnees, les equations des 
10 deux proprietes du revetement, a savoir respecti vement le 
poids de couche pdc et le nombre de cristaux nc par unite 
de surface, s'ecrivent : 

Y pdc = 1,50 + 0.16 X 3 - 0,13 0.079 Xg- 0.070 X^Xg (I) 

15 Y nc = 119 - 39 X 3 -17 X 2 X 3 -14 X^ + 10 X £ X 4 (II) 

EXEMPLE 2 - Determination de 1* equation pr evisionnelle 

representative d'une propriete d'une surface 
metallique traitee dans un bain de conversion 
chimique . 

20 La propriete prise en consideration est toujours 

le poids de couche. 

Pour la determination de 1* equation correspondan- 

te . on peut prendre en compte non seulement des consti- 

tuants chimiques mais aussi des facteurs physiques tels 
25 que la temperature, la Vitesse d' agitation et autres, ou 

bieh des parametres discrets tels que la nature ou l'etat 

de surface d'un substrat. 

Dans le cas present, on a determine 1' equation 

Ypdc a partir de sept variables. 

30 Six de ces variables, respectivement X 1 . X g , X 3 , 

X . X et X , concernent le bain de conversion chimique, a 
4 5 7 

savoir : 

- X 1 : concentration en ions zinc 

- X« : concentration en ions chlorate 

2 

35 - X : concentration en ions phosphate 

- X^ : taux d ' accelerateur (nitrite de sodium) 
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- Xg : temperature cfu bain 

- X ? : valeur de l'acidite libre 

la septieme, X g . representant la concentration du bain de 
conditionnement de surface relatif a l'etape qui precede 
5 Xe traitement de conversion chimique proprement dit. 

En operant comme a l'exemple 1 quant au mode ope- 
ratoire, X' equation Y pdc obtenue s'ecrit : 

Y pdc = 2,15 - 0,19 X_ + 0,16 X. - 0,13 X 0 * 0,13 X^ + 

1Q 0,11 X 3 - 0.10 X 4 + 0,07 X ? (III). 

EXEMPLE 3 - ExempXes pratiques d ' application . 

tin revetement de conversion chimique teX que ceXui 
obtenu avec X* agent de conversion commercialise par la 
Societe Demanderesse sous la marque " GRANODINE 908" doit 
15 realiser 1' adequation suivante : 

son epaisseur ou poids de couche doit etre suf- 
fisamment faible, de l'ordre de 2 a 2,5 g/m , pour permet- 
tre une bonne adherence de la peinture, 

sa structure cristaXXine doit etre Xa pXus dense 
20 et Xa plus reguliere possibXe pour minimiser Xa porosite 
et pour ameliorer Xa protection qu'il con-fere contre Xa 
corrosion . 

C'est cette adequation "adherence peinture-densite 
de structure cristaXXine" qu'iX est propose de resoudre a 
25 l'aide de X* equation indiquee a X'exempXe 2. 

De fa?or> generale. Xes parametres d'un teX bain 
varient suivant Xes fourchettes ci-dessous : 
0,a g/X < Zn ++ < 1,4 g/X 
15 g/X < P0 4 < 23 g/X 
30 0.3 g/X < CX0 3 < 0,9 g/X 

0 . 5 g/X < Ni < 0 , 8 g/X 
0.9 mX NaOH Hf 1 Q < AcL < 1 . 3 ml NaOH N/10 

5Q*C < temperature du bain de conversion < S0*C 
1 , 0 ml < taux d'accelerateur < 2 , 0 ml 
35 T,0 g/1 < concentration du bain de conditionnement 

de surface "FIXODINE 5" < 3,0 g/1. 
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On mesure tous les parametres X^X^ a des instants 
successifs, par exemple a intervalles d'une demi-heure. 
Par exemple : 

a ) Mesure a l'instant : 
Concentration du bain de conditionnemen t de surface: 2 g/1 
Parametres du bain de conversion chimique : 

Zn ++ : 0,8 g/1 

C10 3 - : 0.6 g/1 
P0, 3 ~ : 19 g/1 

taux d ' a cceler ateur : 1,5 ml 

temperature : 55*C 

acidite libre : 0.9 ml. 

On introduit cette serie de valeurs dans Inequation 

(III) Y pdc, apres les avoir centrees et reduites : 

Y pdc =: 2,15 - 0.19 X 6 + 0.16 - 0,13 X g ♦ 0.13 X g + 

0,11 X 3 - 0,10 X A + 0,07 X ? (IV) 

On obtient alors, pour 1' instant t : 

2 

Y pdc = 2.06 g/m . 

Cette valeur de poids de couche qu'on obtiendra se 

2 

situe bien dans la fourchette comprise entre 2 et 2,5 g/m 
qui a ete preconisee precedemment pour l'obtention d'une 
bonne adherence de la peinture. 

II n'y a done pas lieu d* intervener sur les para- 
metres. 

b ) HesMye a \ \ jngtant; t 2 : 
Concentration du bain de conditionnement de surface: 1 g/1 
Parametres du bain de conversion chimique : 

Zn ++ : 1.4 g/1 (X^ ) 

C10 3 - : 0,3 g/1 (X 2 ) 
P0 4 3 " : 19 g/1 <X 3 ) 

taux d ' accelerateur : 1,0 ml 

temperature : 55 *C 

acidite libre : 1,0 ml. 

En introduisant cette serie de valeurs dans l'equa- 

2 

tion (IV), on trouve Y pdc - 2.77 g/m . 

Cette valeur calculee est hors de la -fourchette 
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acceptable. Elle traduit une derive d'un ou de plusieurs 
parametres qu'il convient de corriger pour l'ajustement le 
plus economique. 

Pour cela , on utilisera les projections de l'hyper- 
5 surface a 8 dimensions representant l'equation (III) sur 
les plans corr espondant chacun a deux parametres. Par 
exemple, la -figure S montre, pour des valeurs constantes 
de X 2> X 3 . X^ , X 5 et X ? correspondant a la mesure a 1' ins- 
tant tgi les courbes d'isoreponse, ou courbes de niveau, 
10 dans le .plan X , X , oii : 

I o 

X^-est la teneur en ions Zn* + du bain de conversion. 
X 6 est la teneur en FIXODINE 5 du bain de condition- 

nement de surface. 
Chaque courbe de niveau correspond a une meme va- 
*5 leur du poids de couche Y,(indique en g/m 2 >. Les valeurs 
des autres parametres sent r 



*2 
X 3 
X 4 

O 

X 7 



25 



35 



0.3 g/1 de C10 3 
t9 g/1 de po 4 3 " 



t ml d " accelerateur 
55*C 

1.0 ml d'acidite libre. 
La calculatrice va alors intervenir pour determi- 
ner, a partir des valeurs de X, les variations optimales a 
provoquer dans le bain de conversion chimique et/ou le 
bain de conditionnement pour ramener Y dans la plage ac- 
ceptable, par exemple a une valeur proche de 2 g/m . 

Pour cela. a chaque scrutation. toutes les valeurs 
X sont introduites en memoirs de travail, soit automati- 
quement. soit manuellement a la fin d'une analyse d'echan- 
30 tillon du bain. La calculatrice cherche alors. par la 
methode du gradient. les variations optimales de chacun 
des parametres conduisant a l'obtention de la valeur de 
poids de couche voulue. 

On rappelle que la methode du gradient permet 
d'effectuer des displacements sur 1 " hypersurf ace en chaque 
point perpendiculairement aux courbes isoreponses. le gra- 
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dient etant la direction qui permet de se deplacer le plus 
rapidement vers un optimum. 

On appelle gradient d * une fonction u = f <x,y,z)V 
le vecteur dont les projections sur les axes de coordon- 
5 nees sont les derivees partielles correspondantes de la 
fonction don nee . 

grad u= diii*ilM3 + «Ctuk 
dx dy dz 

i, j, k representant les vecteurs unitaires du systeme de 
reference choisi. 

Le gradient d'une fonction de trois variables est 
en chaque point dirige suivant la normale a la surface 
passant par ce point. La direction du gradient d'une fonc- 
tion en un point donne est la direction suivant laquelle 
la fonction a la plus grande vitesse de croissance en ce 
point, c'est-a-dire pour 1 = grad u la derivee itu atteint 
sa plus grande valeur qui est egale a : dl 

l dx J ^y 1 l dz } 

Oe meme on definit le gradient d'une fonction a n 
variables u = f ( x , y , z , , n ) 

rad u= j4lAi+ilM:S+JiUk + ... + dLun 
dx dy dz dn 

Celui-ci est en chaque point dirige suivant la normale a 
X* hypersurf ace passant par ce point. 

Apres traitement, la calculatrice fournit les 
valeur s optima les 
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x i 

30 X 2 

X 3 
X * 
X 5 
X 6 

35 X 7 



1 , 09 g/1 de Zn + * 

0,56 g/1 de C10 3 ~ 
16 g/1 de PO, 3 ~ 

1 . 33 ml 
51 *C 

2.25 g/1 

0,91 ml 

2 



correspondant a Y - 1,94 g/m (choisie inferieure a l'op- 
timum etant donne la poursuite de Involution). 



24 



0196259 



On supposera maintenant que 1 ' a jus tement doit se 
faire par 

- apport de FIXOOINE 5 dans le bain de conditionnement, 
sous -Forme d'un premelange aqueux contenant 3 g/1 de 
5 constituants actifs, 

apport en produit phosphatant realise a partir de deux 
solutions aqueuses concentrees : 

• l'une A permettant la regeneration du bain de con- 
version en f onctionnement normal, contenant des 
10 ions Zn ++ , P0, 3 ~. Ni ++ , par exemple dans la formu- 

4 

lation suivante {% en poids) 
3.0 < Zn + + < 6.5 

30 < H 3 P0 4 < * 5 
0.1 < Ni** < 1,9 
15 dont l'acidite libre est comprise entre 3 et 11 ml 

de NaOH N (prise d'essai (p.e.)r 1 ml) et l'acidite 
totale est comprise entre 9,8 et 13,4 ml de NaOH N 
(prise d'essai : 1 ml), 
. 1' autre 8 permettant de corriger les augmentations 
ZQ de l'acidite libre et de la concentration en zinc 

du bain et contenant des ions 90. , Ni +f , par 

. 4 

exemple. dans la -formulation suivante : 
30 < NaH 2 P0 4 < 40 
0,1 < Ni* + < 1 
25 4 , 2 < pH < 3,8 

3 < AcT < 5 ml NaOH N (p.e. 1 ml) 
apport en accelerateur s dans le bain de conversion chi- 
mique realise separement a partir de solutions aqueuses 
concentrees : 

30 . l'une C contenant des ions nitrites dans une -Formu- 

lation comprise, par exemple, entre 16 et 20 %. 
. 1" autre D contenant des ions chlorates dans une 
formulation comprise, par exemple, entre 20 et 30Z. 
L'ajustage du bain de conditionnement de surface 
35 peut s'effectuer de fa^on tres simple. Le calculateur 10 
envoie a la carte de commande correspondante une valeur 
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representative de 1 • augmentation de teneurs a obtenir, 
c-est-a-dire 2,25-1= 1, 2 5 g/l. La carte comports un micro- 
processeur ou simplement une memoire morte de codage qui 
determine, a partir de cette valeur et du volume du bain, 
le volume de premelange a ajouter et ouvre une electro- 
vanne a debit constant d'amenee de premelange pendant la 
duree necessaire. 

L'ajustage du bain de conversion est un peu plus 
complexe. La calculatrice determine les apports de solu- 
tions A. B. C et D requis pour amener le bain a la compo- 
sition requise et envoie les informations correspondantes 
aux cartes de commande 56. Chacune de ces cartes comporte 
un circuit de determination de la duree d ' apport necessai- 
re a partir d'un reservoir correspondent 58 et un circuit 
de puissance qui ouvre 1 ' electrovanne appropriee 60 pen- 
dant le temps requis. 

Une autre carte 62. adressable comme les cartes 56. 
repoit I* information de temperature et commande un circuit 
de commutation 64 reglant la puissance electrique appli- 
quee par une source 66 a une resistance chauffante 68. 
dans le sens de la diminution dans le cas considers. 



* 

* * 



II est egalement possible de reguler un bain du 
genre en question en effectuant un compromis entre plu- 
sieurs proprietes. 

Ainsi. on peut utiliser simultanement les equa- 
tions previsionnelles Y pdc (I) et Y nc (II) de l'exemple 
1. Le systeme permet alors de faire evoluer les parametres 
simultanement pour obtenir le meilleur compromis entre le 
poids de couche et la densite de structure. 

Comme il va de soi et comme xl resulte d'ailleurs 
deja de ce qui precede. 1' invention ne se limite nullement 
a ceux de ses modes d ' application et de realisation qui 
ont ete plus specialement envisages ; elle en embrasse. au 
contraire, toutes les variantes. 
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REVENPICATIQNS 
1 . Procede propre a permettre : 

d'une part, de connaitre. a partir de la mesure 
de divers parametres evolutifs d'un milieu chimique evo- 
5 lutif , les particularites du resultat auquel conduit le 
milieu en question a 1" instant choisi, 

d'autre part. au cas ou la valeur trouvee par 
cette mesure montre que le resultat en question ne corres- 
pond pas a celui souhaite, de determiner tout d'abord un 
10 ensemble de modifications des variables ou parametres du 
milieu, propre a faire en sorte que le resultat auquel 
conduira le milieu ainsi modifie corresponde a celui sou- 
haite, puis apporter cet ensemble de modifications au 
milieu chimique evolutif en question, par action sur les 
15 parametres concernes r 

ce procede etant caracterise par le -Fait 

que, pour au moins l'un des resultats auxquels 
conduit ce milieu, on etablit une equation dite prevision- 
nelle du type 

20 Y = + b lXl ♦ b 2 X 2 ... + b k X k ♦ b n X n * b 12Xl X 2 .... 

+ b k-1.k X k-1 X k b k k X k 2 b nn X n 2 

dans laquelle 

- Y represente le resultat auquel conduit le milieu, 

- X, . X . . - . X , ...X representent les variables ou 
zd 1 z k n 

parametres du milieu chimique evolutif dont depend le 
resultat Y et 

- b« r b„ . ..b f b sont des coefficients par 

0 1 k nn 

rapport; auxquels s'exprime la linearite de 1' equation, la 
valeur de ces coefficients etant determinee par des cal- 
culs reposant sur des methodes de regression en elles- 
memes connues, 

que I'on mesure de maniere. repetitive chacun 
desdits variables ou parametres par autant de sondes ou 
capteurs et 

que. d'une part, a un instant t donne , c'est-a- 
dire a 1' issue d'une des mesures repetitives , on calcule 
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35 



27 



0196259 



la valeur de chacune des equations Y en utilisant les va- 
leurs mesurees pour les different* variables a l'instant t 

que, d'autre part, even t uellement et au cas ou 
la valeur trouvee pour lune ou l'autre des equations Y ne 
correspond pas au resultat souhaite, on determine par le 
calcul un ensemble de variations ou modifications a impar- 
tir aux variables ou parametres, cet ensemble de varia- 
tions dormant pour Y la valeur souhaitee et constituant, 
parmi toutes les solutions possibles de I'equation previ- 
sionnelle Y, celle qui est la plus economique et la plus 
facile a realiser. et 

que, d'autre part encore, on agit individuelle- 
ment sur autant d'organes de reglage que de parametres, 
fournissant a chacun de ces organes de reglage un signal 
propre a amener la variable ou parametre qui lui est asso- 
cie a la valeur corr espondan t a la susdite solution de 
1 ' equation pr e vis ionnelle Y la plus economique et la plus 
facile a obtenir. 

2. Procede selon la r evendica tion 1, caracterise 
par le fait que l'on effectue 1* ensemble des calculs a 
I'aide d'une calculatrice numerique que l'on munit de 
bases de donnees appropriees. 

3. Procede selon l'une des revendications 1 et 2 , 
caracterise pax le fait que l'on calcule les modifications 
a impartir aux variables ou parametres par le processus 
represents sur les figures 4a-4c. 

4. Installation permettant de controler un milieu 
chimique evolutif et de le reguler, pour parvenir a une 
valeur optimale d'une fonction predeterminee de parametres 
physiques ou chimiques represent a tifs du milieu, la fonc- 
tion Y etant reliee aux valeurs des parametres par une 
relation de la forme : 

Y ■ b 0 * b 1 X 1 + b 2 X 2 ••• * \\ + b „ X n * > 12 V 2 ' 

* b k-1. k X k -1 X k b kk X k 2 b nn X n 2 

dans laquelle 

- Y represente le resultat auquel conduit le milieu, 
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- Xj t . . . • r . . -X^ representent les variables ou 

parametres du milieu chimique evolutif dont depend le 
resultat Y et 

" b o ' b 1 "" b k b nn sont des coe ^' fi cients par 
5 rapport auxquels s ' exprime la linearite de 1" equation, la 
valeur de ces coefficients etant determinee par des cal- 
culs reposant sur des methodes de regression en elles- 
memes connues, 

ladite installation comprenant : 
10 - des capteurs destines a fournir la valeur ins- 

tantanee de chacun desdits parametres. 

un organe de scrutation permettant de relever 
periodiquement la valeur de 1* ensemble des parametres, 

un organe de calcul permettant de determiner la 
15 modification a apporter auxdits parametres pour donner une 
valeur optimale a la fonction Y, et 

- des moyens de reglage desdits parametres comman- 
des par ledit organe de calcul pour amener la valeur des- 
dits parametres a 1" optimum. 

20 



25 



30 



35 





FIG.2. 
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DEBUT 




Introduire N : nombre de facteurs experimentaux 
X(i) 



Introduire I, J : indices des facteurs X(i) variables 



I 



K = 0 



K = K + 1 



non 




non 



Introduire la valeur affectee 
au parametre X(K) 



20 

7_ 



Introduire l'equation de la fonction-reponse Y : 
N N N 

Y = b Q + I bi X (i) +ZL b n X(i) X(j) + . 
i=1 i=1 j=1 1J 



I 



Introduire les valeurs limites des 
facteurs variables X(I) x(J) : 



X(I) 
X(I) 



MIN 
MAX 



X(J) MIN 
X(J) MAX 



-22 



FIG .3a 
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i 

Introduire un facteur d'echelle pour X(I):X(I) L 

pour X(J):X(J) 

X(i) : X(I) MAX - X ^)MIN 



X(J) L : X(J) MAX - X ( J )MIN 



(P, Q) dimensions de la grille 

1 
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x(D - x(i) MIN 
*(J) = x(j) MIN 



u = 


» 0 


V = 


• 0 



X(I) 
X(J) 



x(D MIN + u = x(i) L 

X(J) MIN + V = X < J )|_ 



Y(U,V) = Y 




V = V + 1 



U = U+ 1 



FIG.3b. 
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j 
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Introduire : 
NB : nombre d' isoreponses 
Yj : valeur initiale de Y 
Y L : pas de variation de Y 
DY : ecart sur Y 



30 



K = 


0 


U = 


P 


V = 


Q 



L V = Q - 1 




U = P - 1 



K = K + 1 



non 




memorisation de la 
valeur Y(U,V) dans 
la grille. 



FIG.3c 
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Introduire N : nombre de facteurs exp§rimentaux X 

i — — 



Introduire equation de la fonction-reponse Y 
N N H 

Y = b o4 < b i X (i) + hi X ^J> 

i=1 i=i j=i *J 

i<j 



Introduire les derivees partielles D[x(i)] de Y 
par rapport a chacun des paramdtres X(i) 

DCX(i)] = 

OA J 



-34 
-36 



Introduire : 

R : rayon du domaine de validity de la fonction Y 
L : facteur multiplicatif tel que 

M(I+1) = H(i) + L * grad Y(I) 
DL : I i mite infer ieure de L 
DY : limite inferieure de l'accroissement de Y 

a partir duquel 1 'optimisation doit s'arreter 



I 



Introduire les coordonnees du point initial [valeur mesuree des X(i)] 

1 



Calcul de Y et de ses derivees au point initial I 

I 

I Introduire JT : nombre maximum d 1 iterations j 



Y(o) = Y 








J = 1 





F\GAo 
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XP (J, 0) = X(J) , 










X(J) = X(J) + L * D[Y(J)] 
XP (J, 1) = X(J) 








non 



J = J + 1 
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Ca lcul de^Y (IT = 1) et de ses d§riv6es partielles 




J 



X(J) = XP(J, 1) j 




FIG.4b. 



envoi des valeurs de 
consigne appropri§es 
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j = t 



T = T + X(J) 2 




non 



J = J+1 - 1 



oui 
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IT = IT+1 



J = 1 



X(J) = X(J) + L * D[X(J)] 
XP(J, IT) = X(J) 




non 



J = J+1 



Calcul de Y et de ses deriv§es 
partielles 

Y(IT) = Y 



T = 0 




FIG 4c 
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Concentration 
en FIXDDINE 5 

(g/0 
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FIG.6 




i / Concentration 



56-c: 




FIG.7. 
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